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Abstract 

Lattice studies of monopole condensation in QCD are based on the construction of 
a disorder parameter, a creation operator of monopoles which is written in terms of 
the gauge fields. This procedure is expected to work for any system which presents 
duahty. We check it on the Ising model in 2d, which is exactly solvable. The output 
is an amusing exercise in statistical mechanics. 
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1 Introduction 

Duality [1,2] is a property of many systems in field theory and in statistical 
mechanics. Their common feature is that, when viewed as Euclidean field theo- 
ries, they have spatial [(d— l)-dimensional] configurations carrying a conserved 
topological charge. The high temperature (strong coupling) phase is charac- 
terized by condensation of these topological structures. In a dual description 
these structures are described by local fields, and their v.e.v. acts as a dual 
order parameter, or disorder parameter. 

In a few models the transformation to the dual is known analytically (the 
Ising model, the U(l) model with Villain action [3], the XY model). In 
other theories duality is expected to be at work, but its presence has only 
been tested numerically [4-9] . The idea is to guess the topological symmetries 
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classifying the extended structures, construct a disorder parameter as the cre- 
ation operator of a structure with non-trivial topology [10], and then measure 
it by numerical simulations [4-9]. 

The approach [10,4-9] relies on symmetry, and is insensitive to the exact choice 
of the action, or to irrelevant terms, contrary to the explicit transformation 
to the dual, which is only possible for specific forms of the action (the Villain 
action in the XY model and in the U(l) model in (3 + l)d). 

A test of the construction is the determination of the critical indices at the 
phase transition by finite size scaling analysis. This has already been done 
for QCD [9], where the natural topological structures in 3d are monopoles. 
A clear evidence has emerged of dual superconductivity as the mechanism of 
colour confinement. 

In (2 + l)-dimensional models the topological structures are vortices (XY 
model) or 2d 0(3) instantons (Heisenberg model). The XY model in its Villain 
form was already known to present duality and condensation of vortices in the 
strong coupling phase. Numerical studies by the symmetry approach confirmed 
that. The Heisenberg model was first viewed as a dual system in Ref. [8]. 

In this paper we test the symmetry approach on the simplest and prototype 
model in (1 + l)d presenting duahty: the Ising model. We show as a theo- 
rem that our construction of the disorder parameter is identical to that of 
Ref. [2]. We give an explicit construction of the dual variables, and we check 
our numerical approach vs the exact solution of the model. 



2 Disorder parameter 



The action of the 2d Ising model is 

S[a] = — J ^ a{n)a{n + (i) , 



where n = (no, rii) is the coordinate of the lattice site and /t is the unit vector 
in the direction — 0,1. cr{n) is the field variable, with values a{n) — ±1. 



The partition function is 



Z — exp 

{a{n)} 



(2) 



The model is exactly solvable [11] and presents a second order phase transition 
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at the critical temperature 



Tc = . (3) 

ln(v^+l) 



The order parameter is the magnetization 

M= lini ^EM^))- (4) 

K— >oo V 

n 

At low temperatures T < Tc, (a) ^ and the symmetry of the action under 
the transformation 

-a (5) 

is spontaneously broken. For T > Tc, {a) = (in the infinite volume limit 
V^oo). 



The critical exponent known as /3, which governs the behaviour of (a) near 

(a) oc (Te - T)^ , (6) 



T 



has the value 



(7) 



The correlation function of the field a 

r{%j) = {a{t)a{j)) (8) 
near the phase transition is approximately rotation invariant, only depends 



on the distance d — yit — j)"^ and behaves as 

^~exp(^-^j+(a)^ (9) 

For T > Tc, (cr) = 0, and the correlation length ^(T) diverges at T = Tc as 

e = a(T-Tc)-V (10) 
The critical index u has the value u — 1. 

A dual description of the system can be given in terms of a dual field n, 
defined on the dual lattice [2]. The dual lattice associates to each plaquette 
of the original lattice a site, which can be identified as its centre, and to each 
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link a dual link, which is orthogonal to it and cuts it. The variable // is defined 
throughout its correlation functions 



(a^(7)/^(J)) = §11 ■ (11) 

S differs from S by the introduction of a magnetic dislocation, i.e. by changing 
signs to J on all the links crossed by a line joining ?and Jon the dual lattice: 
the result is independent of the choice of the line. 

The field fj, obtained in this way is again valued /i — ±1 and the partition 
function Z[a, T] obeys the equality 

Z[a,T] = Z[f,,T*] , (12) 

with 

2 1 

sinh — = — - — 7y- , (13) 
T sinh^ ' ^ ^ 

i.e. T* ^ Tc if T ^ T^. 

The dual description of the system maps the disordered phase (a) — onto the 
ordered phase of the dual. For T > Tc, (/i) ^ 0, (a) = 0. The order parameter 
of the disordered phase is called a disorder parameter. Near the critical point 

(/.)«(T-Te)\ (14) 

which defines the critical exponent S associated to the disorder parameter. 
Then, self-duality of the model implies that 

S = p^l. (15) 



We shall approach the problem in a slightly different form, emphasizing the 
symmetry aspects of the disordered phase. We shall then prove that our dis- 
order parameter is equal to that of Ref. [2]. A similar procedure was used 
in Ref. [6] for the U(l) gauge theory, as an alternative procedure to that of 
Ref. [3]. 

The Ising model can be viewed as a 1 -|- 1 dimensional Euclidean field theory, 
with action 

^ M=0,1 n 
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with A^(7(ri) = a{n + ji) — a{n). The current 

= ^e^i/A^cT (17) 

is identically conserved, 

^,J, = . (18) 
The corresponding constant of motion is 

Q = J2Mno:ni) (19) 
m 

or 

1 1 
Q = Ai(T(no,ni) = - [(T(no,+oo) - (T(no, -oo)] . (20) 



Q is a topological charge, which classifies 1-dimensional spatial configura- 
tions by their boundary conditions. A kink corresponds to a{no, +oo) — 1, 
(j(no, — oo) = — 1 or Q = 1, an antikink to Q = — 1- 

Our guess is that the conservation law Eq. (18) is spontaneously broken in the 
high-temperature phase, by condensation of kinks. As a disorder parameter 
we shall choose the vacuum expectation value of an operator carrying nonzero 
charge Q (Eq. (19)). For that operator we choose the creation operator of a 
kink. 



The creation operator of a kink at site rii and time no is defined as 
n{no,ni) = exp 



-2^ ^ a{no,n)a{no + l,n) 

n<ni 



(21) 



The correlator 

P(no) = (/i(no,ni)//(0,ni)) (22) 

is accordingly defined as 

1 / S'\ Z 

V{no) = - ^ exp (^-- j /i(no, ni)/i(0, m) = - . (23) 

Z is obtained from Z by reversing the sign of the temporal links Ja{Q, n)a{l, n) 
and Ja{no, n)a{no + l,n), n < rii, in the action. 

It is trivial to check that 'D^Uq) Eq. (22) really corresponds to the propagator 
of a kink at the spatial site ni from time t = to time t = Uq. Indeed in com- 
puting Z we can change variables in the sum from a{l,n) to — cr(l, n), n < rii. 
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This change brings the temporal hnks 0-1 to the original form they appeared 
in Z , but changes sign to the spatial link o{\. nx)o{\, n\ + 1) and to the tem- 
poral links (t(1, n)(T(2, n), n < n\. A kink has been added to the configuration 
at time Uq — 1. We can now change variables by sending cr(2, n) — > — cr(2, n). 
The result will be again the appearance of a kink at no = 2, and a change 
of sign of links (j(2, n)(T(3, n), n < ni. The construction can be repeated and 
finally the change of sign of the temporal link at no will be reabsorbed by the 
kink sitting at Uq in Z. This completes the proof. 

The net effect is a change of the spatial hnks {ni,ni + 1) at all times in the 
interval (0, no] . This is exactly the definition of the dual correlator given in 
Ref. [2]. So our (/i) is equal to the disorder parameter of Ref. [2] and provides 
an explicit construction of the operator in terms of the original fields a. The 
construction is exactly the same which was given in Refs. [6-8], respectively 
for the U(l), the 3d XY model and the Heisenberg model. 

Some remarks are necessary at this point. First, in the absence of the second 
kink at no, the construction sketched above includes a change of sign of the 
temporal link at the border of the lattice and the periodic boundary conditions 
are made antiperiodic by the change of variables. Second, for finite lattices 
the correlator Eq. (22) is a four point function. As no goes large, by cluster 
property 

V{no) {ii)' . (24) 

Finally, we stress that the result of our construction is mathematically iden- 
tical to that of Ref. [2]. (//) 7^ 0, T > Tc is a theorem. However for more 
complicated choices of the action, belonging to the same universality class, 
our construction, which only relics on symmetry, can be more practical for 
numerical approaches. This is similar to the situation for U(l), where the con- 
struction of the dual is feasible for the Villain form of the action, but the 
approach based on symmetry allows treatment of Wilson and other forms of 
the action. 

As a check of our arguments we shall compute the disorder parameter numer- 
ically, and use it to determine the critical indices of the model. 



3 Numerical simulations 



Let us try to extract information about the phase transition from the numerical 
study of the disorder parameter. To this purpose, at fixed ni in a finite lattice 
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with periodic boundary conditions, we study 

V{no) = I (25) 

as a function of /3 = 1/T by standard Monte Carlo methods. 

By the cluster property, at large time separation it is expected that 

V{no) ^ , (26) 

being (yu) the disorder parameter. Hence in principle we can extract (/i) by a 
fit of 'P(no) data at large Uq. (/j) is expected to have the opposite behaviour 
to the magnetization M: in the thermodynamic limit (/i) is rigorously zero for 
P > Pc: different from zero for /3 < /3c, and for /3 < /9c, it approaches zero with 
the power law 

(/.) ^ t' , (27) 
where t is the reduced temperature {T — Tc)/Tc. 

A direct determination of 'D{no) by Monte Carlo techniques requires an un- 
bearable computational cost. This is a consequence of the form of V^uq) in 
terms of the order variables [4]. In fact ©(no) is the ratio of two partition func- 
tions and the determination of the partition function of a system by numerical 
simulations is very difficult: Z is the exponential of a quantity fiuctuating like 
the square root of the volume, so we expect for it fluctuations of order . Be- 
cause of these strong fluctuations the direct determination of ^'(no) by Monte 
Carlo techniques requires an amount of statistics which can not be obtained 
in realistic simulations. 



This difficulty can be overcome by defining another quantity whose fluctua- 
tions are proportional to the square root of the volume and which gives in- 
formation about the behaviour of (//) . To this goal we introduce the function 

[4] 

p=^logD(no). (28) 

In the limit no — > oo 

p^4^1og(/.), (29) 

which implies 

{l^) = exp(^j%{P')dp') , (30) 
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Fig. 1. /9 as a function of (3 for different lattice sizes. Continuous lines refer to low 
temperature calculations. 

and the behaviour of (n) can be reconstructed from the behaviour of p in a 
simple way. 



By using Eq. (25) we get 



P = {S)s - {S) 



s ) 



(31) 



where (•S')^ is the average of the action of the Ising model, Eq. (1), taken on 
configurations generated with that action, and (S) § is the average of the kink 
action S taken on configurations generated with this modified action. The 
quantity p will be much easier to measure in a numerical simulation, since it 
is expected to fluctuate as W. 



We have investigated lattices of size L x 2L, with the spatial extension L 
ranging from 80 to 200. In our simulations we chose Ui = L/2 and uq = L, and 
periodic boundary conditions were imposed. We simulated both the standard 
and the kink action of the model by means of a heat-bath algorithm. Errors 
were estimated by applying the jack-knife method to bunched data. Far from 
criticality we collected about 100000 measurements for each /3 by sampling 
each two sweeps; near the critical couphng the statistics has been increased 
by a factor of 25. 

From the behaviour of (/x) in the thermodynamic limit, we expect that p is 
nearly constant at low j3 (due to the slow variation of (p) in this region), 
has a sharp negative peak in the critical region (corresponding to the sudden 
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Fig. 2. Low P data for p. 

decline of (/i)) and approaches — oo at weak coupling (as a consequence of the 
restoring of the dual symmetry). The extrapolation of numerical results to the 
thermodynamic limit should give this behaviour for p. A plot of p data for 
the lattices we have investigated is reported in Fig. 1. We shall show that the 
obtained shape of p is compatible with the expected behaviour in the infinite 
volume limit. 

By simple algebraic manipulation, 

V{no) = (e-^(^-^))5 , (32) 

where the subscript S means that the average has to be taken over configura- 
tions weighted with the standard Ising Boltzmann factor. Only modified links 
(i.e. links whose sign has been changed) contribute to — S", so we get 

P(no) = (e-2^i:.'-''^'=)s, (33) 



where at and ak are spins connected by a modified link I' and the sum is 
performed over all modified links. 

In the limit [3 ^ oo the system is completely ordered, so all links are equal 
to 1 and the sum in Eq. (33) simply reduces to the total number of modified 
links, i.e. 

V{no) e-2/^^ . (34) 
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Fig. 3. Rescaled plot of p data. 



Indeed at high j3, (/i) goes to zero in the infinite volume limit, as required 
from the dual symmetry restoration. 



From Eq. (34) we get that 



-2L 



(35) 



at high (3. In Fig. 1 continuous lines indicate the expected results from the 
high (3 calculation. The agreement with numerical results is evident. We can 
also consider the opposite implication: if Eq. (35) holds in the weak coupling, 
in the thermodynamic limit (fi) = in this region. 

On the other hand, p is compatible with zero for a wide range of /5's below 
(3c, independently of the volume (see Fig. 2). This implies that (n) is nearly 
constant in this region and that feature seems to extrapolate to the thermo- 
dynamic limit. Since at high temperature 



(/i) = 1 , (36) 

our data for p imply that (/i) in the thermodynamic limit is different from 
zero in a wide range of /3's below /5c. Hence the behaviour of p we obtained is 
compatible with dual symmetry breaking for (3 < (3c- 



In the critical region the approach to the thermodynamic limit is governed by 
finite size scaling theory, which states that at finite size L in the critical region 
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0.438 


92.3(1) 


80,120,160,200 


0.439 


138.9(2) 


120,160,200 


0.44 


238.1(4) 





Table 1 

Correlation length in the 2d Ising model measured in a 400 x 400 lattice at the (3 
values used in Fig. 3. The third column indicates the lattice sizes which are expected 
to show scaling behaviour (see text). 

for T > Tc, {jj) is given by 

W = t'f{s) , (37) 
with / an unknown function of the scaling variable 

s = tL"" . (38) 

Prom Eq. (37) it follows 

_ A5 4 d , 

Prescal - - ^ ^^^^ ^J^S) , ^6^)) 

i.e. the rescaled variable (p/3^)/(L^/'^/3c) is a function of the scaling variable s. 

A plot of rescaled data is shown in Fig. 3. According to Eq. (39), all points 
should fall on one single curve. We observe however deviations at small and 
large s. It is important to understand the range of validity of Eq. (39). On the 
one hand, in order for Eq. (26) to be valid, the antikink should be placed at 
a large temporal distance from the kink. This means that 

no^L>i, (40) 

where ^ is the correlation length at a certain temperature. On the other hand, 
we want to be in the scaling region where finite size scaling is valid, so L/^ 
should not take very large values. In order to have an approximate idea, we 
computed the correlation length in a large LxL lattice, L = 400, at the values 
of /3 used to obtain Fig. 3. The results are shown in Table 1. 

The condition Eq. (40) means that, for example for the L = 80 lattice, only the 
points of /9 < 0.438 from the list shown in Table 1 should give a proper scaling 
of p. On the other hand, L should not be much more larger than ^, because 
otherwise we go out from the finite size scaling region. We can tentatively 
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6-20 


0.120(5) 1.42 


0.135(2) 1.21 


0.132(10) 


0.56 


7-20 


0.120(6) 1.53 


0.136(2) 1.22 


0.128(9) 


0.36 


8-20 


0.110(3) 0.52 


0.143(3) 0.48 


0.120(13) 


0.36 



Table 2 

Different fits to obtain the 6 exponent from Eqs. (41)-(43) (see text). 

write L/C, ^ 3, though this value is somewhat arbitrary. In this case this 
choice corresponds to eUminate the final four points in Fig. 3. Condition (40) 
eliminates the five points closer to s = in that figure. Putting the two 
conditions together, the resulting lattices which should show a proper scaling 
are given in the third column of Table 1. These are the middle section, 6th-20th 
points in Fig. 3, which indeed seem to lie on a single curve. 

As a further check, we can fit the critical exponent 6 by means of Eq. (39) 
using as an input u and f3c {'^ = ^, = 0.44068...). For the unknown term 
depending on /(s) we have several possibilities. We can guess that it is roughly 
constant inside the range of the explored s, and do the fit to (fit (a)) 

Prescal ; \^^) 

S 

or consider that when s —>■ 0, both /(s) and its derivative go to a constant, 
and write (fit (b)) 

45 

Prescal = + C , (42) 

S 

where C is a constant term. We can finally use the three-parameter fit (fit (c)) 

_ 46 {s + g) 



The results of the fitsQ are given in Table 2. Eliminating the 6th or the 7th 
point (second and third rows in Table 2) corresponds to taking Eq. (40) as a 
rigourous inequality, L > C^, according to the data of Table 1. Fit (c) gives the 
more stable values with respect to the modification of the number of points 
in the fit, and also those with less by degree of freedom (DF). A mean of 
the values given by the fit (c) (third column in Table 2) gives 6 = 0.127(8), in 
fair agreement with the expected value 6 = 0.125. This fit for points 7-20 is 
shown in Fig. 3. 

^ The fits have been performed by using the MINUIT program. 
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A more complicated fit allowing also the determination of v and /5c in principle 
is also possible, but an appropriate guess of the form of / is needed (see [6- 
8]). This fit is useful when only approximated values of /3c and v are known. 
In our case the solution of the model is known and the simple quality of the 
scahng is itself a good indication of the correctness of the approach. We note 
however that in order to obtain an accurate result for the critical exponents, 
one should not approach very much the critical point, where the correlation 
length diverges, in contrast with ordinary analyses. This can be taking into 
account by monitoring the correlation length, as we showed in this simple case. 



4 Conclusions 

The analysis of the Ising model suggests that the study of p is suitable for 
investigating order-disorder phase transitions from a dual point of view. While 
(/x) is not adequate for numerical simulations, p data are completely reliable 
and allow an unambiguous reconstruction of the shape of (yu) . In addition, from 
a finite size scaling analysis we can get the critical exponent associated to the 
disorder parameter, the critical temperature and the critical exponent of the 
correlation length, though a careful analysis for this observable is convenient 
in order to get an accurate measure for the critical exponents. 

This analysis is the first complete test of the procedure used to demonstrate 
that dual superconductivity of the vacuum is the mechanism of colour con- 
finement. The result of this test confirms the hypotheses that are at the basis 
of the study of Ref . [9]. 
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